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ABSTRACT 

Coronagraphic imaging with the Near Infrared Camera and Multi Object Spectrometer on the Hubble Space 
Telescope reveals a large, ^400 AU (4") radius, circumstellar disk around the Herbig Ae/Be star HD 141569. A 
reflected light image at 1.1 //m shows the disk oriented at a position angle of 356 ± 5° and inclined to our line 
of sight by 51 ± 3°; the intrinsic scattering function of the dust in the disk makes the side inclined toward us, the 
eastern side, brighter. The disk flux density peaks 185 AU (1 ."85) from the star and falls off to both larger and 
smaller radii. A region of depleted material, or a gap, in the disk is centered 250 AU from the star. The dynamical 
effect of one or more planets may be necessary to explain this morphology. 

Subject headings: circumstellar matter — stars: emission line — stars: individual (HD 141569) 



1. INTRODUCTION 

The discovery of disks around main-sequence stars has 
opened possibilities for studying nascent planetary systems. 
Old Herbig Ae/Be and (3 Pic analog stars lie at a transitional age 
between pre- and zero-age main sequence stars. Both classes 
have infrared excesses indicative of the presence of warm cir- 
cumstellar dust as well as luminosities and colors consistent 
with those of young stars. Their circumstellar material, in the 
form of disks, may be composed of both remnant proto-stellar 
and secondary "debris" dust generated by collisions. 

HD 141569 (= IRAS 15473-0346= SAO 140789), classified 
as B9.5Ve by Jaschek & Jaschek (1992), is a Herbig Ae/Be star 
with double peaked Ha emission. Jaschek, Jaschek, & Egret 
(1986) first identified it as a main sequence star with infrared 
excess and an SED peaking at ^60 /im. The disk has an optical 
depth Ti-iaQnm = Ldisk/L* = 8.4 x 10~ 3 (Z uckerman, Forveille, 
& Kastner 1995), which is a few times larger than the two well- 
known infrared excess stars, Pic and HR 4796A. With an Hip- 
parcos measured distance of 99 *g pc (so 1" « 100 AU), the star 
may lie on the out skirts of the dark cloud complex L134/L134N 
( [Sahu et al. 1998 ). Assuming its B-V = 0.095 mag color is in- 
dicative of interstellar reddening, HD 141569 has an absolute 
magnitude My=l-6 mag (van den Ancker et al. 1998). Such 
a low luminosity for its spectral typ e is consistent w ith that of 
other young main sequence A stars (Jura et al. 1998). 

The age of HD 141569 is not well determined. It is located 
approximately on the zero age main sequence, has an infrared 
optical depth similar to that of other young A stars, h as strong 
CO emission ( Zuckerman, Forveille, & Kastner 1995), an d has 
two T-Tauri like companion stars ( Weinberger et al. 1 999 ). To- 
gether, these constrain its age to 10°— 10' yr. 

The 1.1 /im results reported here are part of our guaranteed 
time project to study the environments of nearby stars with the 
Near Infrared Camera and Multi Object Spectrometer (NIC- 
MOS) aboa rd the Hubble Space Telescope (HST) (S chneider 
et al. 1998). HD 141569 was included because of its large 
far-infrared optical depth, CO emission, and extended elliptical 



disk in mid-infrared images by Silverstone et al. (1999). The 
scattered light disk around HD 141569 has also been observed 
with NICMOS at 1.6 /xm by Augereau et al. (1999a). 

2. OBSERVATIONS 

Coronagraphic observations of HD 141569 were made on 27 
September 1998 with NICMOS Camera 2. The F110W filter 
(A CCT!fra /=1.10 fj,m, FWHM=0.59 /im) was used to maximize the 
angular resolution while minimizing scattered light. With the 
target positioned in the 0."3 radius coronagraph, three MULTI- 
ACCUM sequences were obtained, for a total integration time 
of 607.9 sec at each of two spacecraft orientations, rolled 8° 
about the target axis with respect to each other. To minimize 
time-dependent point spread function (PSF) variations, the two 
image sets were obtained within 40 minutes. 

Prior to each set of coronagraphic observations, 0.36 s target 
acquisition frames were obtained with the F171M (A ra „ ra /=1.72 
/im, FWHM=0.07 /im) filter. Contemporaneous lamp flats and 
backgrounds were obtained at F160W for the purpose of locat- 
ing the coronagraphic hole and enabling good flat-fielding near 
the hole. 

3. DATA ANALYSIS 

The MULTIAC CUM data sets w ere processed with the Ni- 
creduce software (MacLeod 1997). Dark frames taken on 20 



September 1998 were used to subtract the dark current and cor- 
rect the detector shading. A special coronagraphic flat field 
was created by modifying a standard high S/N on-orbit Fl 10W 
reference flat field. Within a radius of 12 pixels (0."91) of 
the hole, the pixel sensitivities were measured in the F160W 
hole-finding image, scaled to Fl 10W, and replaced in the refer- 
ence flat. After bias subtraction, linearization and flat fielding, 
known bad pixels were corrected with a distance weighted av- 
erage within a radius of eight pixels, and the three MULTIAC- 
CUM images from each spacecraft orientation were medianed. 
The best available photometric calibration was applied to the 
final Fl 10W images in which 1 ADU s" 1 = 2.031 x 10" 6 Jy and 
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mag = 1775 Jy ( [Rieke 1999| ). 

In the reduced and calibrated NICMOS images, scattered 
light from the star still dominates the flux levels at radii of 
1 — 5", so in order to see the disk, a PSF star must be sub- 
tracted. The same observing strategy and the Fl 10W filter were 
used for 14 other stars in this GTO program, and targets which 
do not show extended disk emission form a library of possible 
PSF stars. The detailed character of the NICMOS PSF changes 
with time pri marily due to thermal variations in the HST opti- 
cal assembly (Kulkarni et al. 1999); hence some PSF star sub- 
tractions produce much lower residuals than others. Each PSF 
image was registered using cubic convolution interpolation to 
both orientations of HD 141569. Only stars brighter than HD 
141569 were considered, and the PSF subtraction residuals are 
dominated by systematic uncertainties. 

Scattered light from two nearby stars outside the field of 
view, see §5.1, pollutes the northwest corner of the image. For 
both roll angles, we used synthetic PSFs generated with Tiny 
Tim ( Krist & Hook 1997 ) and positioned them at the locations 
of the two stars to subtract this excess scatter. A small amount 
of residual flux, however, can still be seen in Figure |jj 

Figure || compares the azimuthally averaged disk flux ob- 
served after subtraction of each of the three best matched PSF 
stars. Two of these PSF stars, t 1 Eri and 49 Cet, were disk can- 
didates in which no disk was detected, and the third, HR 4748, 
was observed specifically as a PSF star. At radii greater than 
l."6 (160 AU), the morphology of the disk is consistent in all 
three subtractions. Subtraction of 49 Cet yielded the smallest 
amplitude residuals, and all further analysis was carried out on 
images using it as the PSF. 

The PSF-subtracted images from each telescope roll were ro- 
tated to a common orientation and averaged. Regions of the im- 
age obscured by the diffraction spikes in one roll were replaced 
where possible by uncorrupted regions in the other roll. 

4. RESULTS 

The PSF subtracted, roll-combined coronagraphic image of 
HD 141569 is shown in Figure Ilk. An elongated disk is clearly 
evident at radii from ~ 1.6 to 4 / '(160 to 400 AU) from the star. 
Elliptical isophotes were fit to the disk image yielding a posi- 
tion angle of 356 ± 5° and ratio of the minor to major axes of 
0.63 ± 0.04. If, as is likely for dust orbiting a star, the disk 
is intrinsically circular, the axial ratio implies an inclination of 
51 ±3° from face-on. This geometry was used to "deproject" 
the disk to a face-on view (see Figure [j]b). 

The azimuthally averaged disk surface brightness, shown in 
Figure ^ as computed from the deprojected image, peaks at 0.3 
mJy arcsec -2 at a radius of 185 AU and falls off to larger and 
smaller radii. Fitting the slope of the surface brightness with a 
power law gives r~ 32 for 190 < r < 250 AU. At 250 AU an an- 
nulus of lower surface brightness, or a "gap" in the disk, can be 
seen, which is apparent in Figure || where the surface brightness 
has been multiplied by 4ot 2 to determine the surface density of 
scatterers (see § 5.2). On either side of this gap, the surface den- 



The corrected total flux density is 8 ± 2 mJy. The Fl 10W mag- 
nitudes of 49 Cet and HD 141569 were extrapolated from their 
H-band magnitudes, and the uncertainty in their flux ratio com- 
pletely dominates the uncertainty in the photometry of the disk. 
A change of ~2% in the subtracted PSF flux, makes a ^20% 
difference in the integrated disk flux density. 

The eastern side of the disk is 1.5±0.2 times brighter than 
the western side. This difference is not an artifact of the PSF 
subtraction and can be explained by dust particles which prefer- 
entially scatter in the forward direction. Such phase functions 
apply even when the particles are a few times larger than the 
wavelength of the scattered light (van de Hulst 1957) and are 
produced in rec ent models of fluffy grains (e.g. W olff, Clayton, 
& Gibson 1998). Sma ll particles probably do exist in the HD 
141569 disk (see §5.2). The brighter side, therefore, is inclined 
toward us, and the ratio of forward to backward scattering im- 
plies an asymmetry parameter, g, of 0. 1 1 in the analytic phase 
function of Henyey & Greenstein (1941). This is consistent 
with the results of Augereau et al. (1999b) on HR 4796A in 
which the upper limit for g was 0.15. 

5. DISCUSSION 
5.1. Disk Morphology 

The surface brightness of the disk has a complicated radial 
morphology which cannot be represented by a single power- 
law. The density of scatterers, assuming the disk is everywhere 
optically thin, is as high at radii up to about 360 AU from the 
star as it is in the region 200 AU from the star. 

Dynamical processes could play a large part in shaping the 
disk. HD 141569 may be the brightest member of a multiple 
star system, since two bright stars (henceforth "B" and "C") 
with H = 8.6 and 9.4 mag lie 7."57 (^750 AU) and 8."93 (^900 
AU) away at respective position angles of 311.5° and 310.0°. 
The high probability that these stars are physical companions 
and assessments of their ages are discussed in a separate paper 
(Weinberger et al. 1999). If the three stars are physically as- 



sity of scatterers does not vary appreciably. At r;>340 AU, the 
surface brightness of the disk falls off more steeply, as r" 5 9 . 

The disk flux density, in a rectangle of area 63 sq. arcsec out- 
side of 0."6 from the star and extending to 5. "3 from the star 
along the major axis of the disk, is 7.2±1.8 mJy. This integra- 
tion excluded pixels in the diffraction spikes; we make a purely 
geometric correction by finding the average flux density at ev- 
ery radius and multiplying it by the number of excluded pixels. 



sociated, the two companions would likely have a significant 
dynamical effect on the disk. If the companion stars are in the 
plane of the disk, the physical separation of "A"-"B" is 990 AU 
and that of "B"-"C" is 190 AU. Then, the ratio of the semi- 
major axes of the wider to the closer pair would be only ~5.2. 
This is not ex pected to be a stable triple system (E ggleton & 
Kiseleva 1995), although such young stars may not yet have 
had time to become unbound. 

If the companions are out of the disk plane, they could be 
much further from the primary, and their dynamical effect could 
be to excite significant velocities in the disk dust perpendicular 
to the plane. This would puff up the disk and possibly flare 
it. Increased collision velocities could promote the breakup of 
larger particles into smaller and could help account for the large 
number of scatterers far out in the disk. 

The dip in surface brightness at r«250 AU which is most 
easily seen along the major axis (Figure [l]), is a striking feature 
of the disk. The full width at half minimum of this gap is 60 
AU, to which the intrinsic size of the NICMOS PSF contributes 
~5%. The depth of this dip in terms of the density of scatterers 
(tw; Figure g) is ~14% of the average value between 185 and 
340 AU. 

No point source is seen within the disk, and to constrain 
the luminosities of possible comp anion bodies, we plan ted PSF 
"stars," generated with Tiny Tim (Krist & Hook 1997), at ran- 
dom locations in the image at one telescope roll angle sub- 
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tracted from the image at the other roll. The limiting magni- 
tude (3cr) of an object which could be detected at each radius 
is plotted in Figure This limit does not apply to the region 
obscured by the diffraction spikes which, at the radius of the 
gap, hide ~17% of the disk. 

At the radius of the gap, the limiting magnitude of 20.3 and 
an assumed age <10 Myr correspond to a planet with mass 
<0.003 Mq (3 Mjup) according to cooling curves by Burrows 
et al. (1997). At this age, the brown dwarf-planet transition oc- 
curs at log(L/L Q )w -3 which corresponds to an Fl 10W magni- 
tude of — 16. Figure |] shows that either a star or brown dwarf 
would be clearly detectable in our image at any radii greater 
than -80 AU. 

If the gap is cleared only by the gravity of a co-orbital planet, 
the planet mass is: M/Mq = C(Aa/a) 3 where Aa and a are 
the widt h and radius of the gap, respectively, and the constant 
C« 0.1 ( pssauer 1993j ). This corresponds to a planetary mass 
of —1.3 Jupiters, which given the detection limits discussed in 



§5.1, could not have been detected. The orbital period for ma- 
terial 250 AU from the 2.3 M Q star is 2600 yr. If we assume 
that it takes >300 orbital periods to clear material in a gap, 
based on models of disks by Bryden et al. (1999), then a gap 
can be opened in ;>8 x 10 5 yr. That the gap is not completely 
clear may mean that processes are moving particles through the 
disk. The Poynting-Robertson timescale for particles J> 1 /im is 

— 4xl0 6 yr, similar to the age of the star. Radiation pressure 
will preferentially clear small grains on a timescale of — 10 4 yr 
and could fill in a gap if such grains are present, or being cre- 
ated, in the disk. 

5.2. Constituent grains 

We can combine the information we infer from the NICMOS 
images with what is known about the disk from mid to far- 
infrared measurements. For reflection from optically thin dust, 
tui = 4-7T0 2 j, where <f) is the angular distance of the scatterers 
from the illuminating star, S is surface brightness, F is the re- 
ceived flux from the star, r is the optical depth of scatterers, 
and u) is the albedo (Figure Bj). At the peak surface brightness, 
0=1. "9, 5=0.3 mJy arcsec"% and ra = 4x 10" 3 . The optical 
depth to visual absorption implied by the far infrared excess is 
8.4 x 10~ 3 ; if this is also the absorption optical depth at 1.1 /xm, 
then the albedo of the scatterers is —0.3. However, about l /$ of 
the luminosity emitted fr om 3.5 - 25 um arises from a region 
within 150 AU of the star dSilverstone et al. 199^). So, the total 



infrared r overestimates the optical depth in the exterior part of 
the disk seen in reflected light and results in an underestimate 
of the albedo. In the region at — 1."9 (190 AU), therefore, the 
albedo is greater than 0.3 and probably more like 0.4. This is 
a substantially hig her albedo than the zodiacal dust in o ur own 
Solar system, 0.1 (Dumont & Levasseur-Regourd 1988), but is 
similar to the albedos of mineral mixtures and ices found in ISO 



spectra of Herbig Ae/Be disks ( jMalfait et al. 1998[ ). 

At r<150 AU, the dust temperature is — 150 K as estimated 
from mid-infrared colors and luminosity (S ilverstone et al. 
1999). If the grains are icy, the decline in 1 . 1 /im surface bright- 
ness inside 185 AU may be due to the sublimation of the ice at 
temperatures —100 K. In this case, the grains emitting in the 
mid to far-infrared are the iceless counterparts of the grains at 
higher radius. 

Gaseous CO was detected by Zuckerman, Forveille, & Kast- 
ner (1995) with a velocity width of 7.6 km s" 1 , corresponding 
to a radius of 86 AU given a disk inclination of 5 1 and stellar 
mass of 2.3M Q . The gas, therefore, is well inside the brightest 
part of the reflected light disk and the gap. In the spectrum of 
Sylvester et al. (1996), a large rise in the continuum shortward 
of 8 /im and a possible peak at 1 1.3 /jm both indicate emission 
from polycyclic aromatic hydrocarbons. The existence of gas 
and small grains in the presence of forces which remove them 
very quickly indicates they are constantly being resupplied. The 
disk seems to be in an active planetesimal building phase. 

6. CONCLUSIONS 

NICMOS images reveal a large circumstellar disk around the 
Herbig Ae star HD 141569. The structure of the disk, including 
a region of depleted material at 250 AU can most easily be un- 
derstood as due to dynamical sculpting by one or more planets 
orbiting within it. If planets exist in this system, they must have 
formed in <10 7 yr and far from the central star compared with 
the planets in our own Solar system. 

We thank David Koerner, Andrea Ghez and Mike Jura for 
helpful conversations. This work is based on observations with 
the NASA/ESA Hubble Space Telescope, obtained at the Space 
Telescope Science Institute, which is operated by the Asso- 
ciation of Universities for Research in Astronomy, Inc. un- 
der NASA contract NAS-26555 and supported by NASA grant 
NAG5-3042 to the NICMOS instrument definition team. 
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FIG. 1 . — Left (a): Image of the reflected light disk around HD 141569 shown in a stretch which is linear with flux density. The eastern side of the disk is brighter 
because of the intrinsic scattering phase function of the dust. Right (b): Deprojected image of the disk using an inclination angle of 51°. Artifacts from small 
instabilities in the PSF subtraction process include the radial "streamers" and the dark "finger" in the inner disk at the 4 o'clock position. These residuals do not 
significantly affect the radial profile or the integrated flux density. The diffraction spikes narrow to higher radius as the overlap between the two sets of observations 
at different orientations increases and allows pixels behind the spikes in one image to be replaced by pixels from the other image. 
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FIG. 2. — Azimuthally averaged surface brightness profile of the disk 
after subtraction of three different PSF stars. The error bars represent only 
the statistical uncertainty in the mean surface brightness at each radius. At 
radii greater than ~160 AU the three profiles are consistent to within these 
uncertainties. The subtractions with r 1 Eri and 49 Cet suggest that the disk 
flux falls off substantially inside of 160 AU. 
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FIG. 3. — The surface density of scatterers plotted as a function of radius 
computed by multiplying the surface brightness by 47rr 2 and normalizing to 
the total stellar flux density. 
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FIG. 4. — Apparent Fl 10W limiting magnitude of an unseen companion 
as a function of distance from the star. The uncertainty in the detection limit 
is the size of the points, or ~10%. 
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